Motility of small nematodes in wet granular media 
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Abstract. - The motility of the worm nematode Caenorhabditis elegans is investigated in shallow, 
wet granular media as a function of particle size dispersity and area density (0). Surprisingly, we 
find that the nematode's propulsion speed is enhanced by the presence of particles in a fluid and is 
nearly independent of area density. The undulation speed, often used to differentiate locomotion 
gaits, is significantly affected by the bulk material properties of wet mono- and polydisperse 
granular media for (f) > 0.55. This difference is characterized by a change in the nematode's 
waveform from swimming to crawling in dense polydisperse media only. This change highlights the 
organism's adaptability to subtle differences in local structure and response between monodisperse 
and polydisperse media. 
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Introduction. — Many live organisms evolve within 
complex fluidic environments such as intestinal fluid, hu- 
man mucus, and soil [T]-[3]. Such organisms often develop 
unique methods of locomotion that allow them to exploit 
the material properties of the media For example, 

in the case of freely swimming spermatozoa in a Newto- 
nian fluid, the fiagellum shows a regular sinusoidal beating 
pattern in both space and time. However, once the organ- 
" ism encounters a complex fluid, like semen, this regular 
beating pattern is replaced by high- amplitude, asymmet- 
ric bending of the fiagellum [6 . This 'hyper- activated' 
sperm is clearly affected by its fluidic environment, which 
in turn, can affect human fertility [T|. 

There has been a renewed interest in swimming at low 
Reynolds numbers p^lTHTT]. where live organisms must 
break time-reversal symmetry [12] in order to achieve net 
motion. The majority of previous work on locomotion, 
however, has been restricted to Newtonian liquids, and few 
investigations are available on the motility of live organ- 
isms in complex fluids such as granular media [2[l4l[5lfT3] . 
Recent investigations [4|E] have indicated that in gran- 
ular media, an organism needs to tune its swimming am- 
plitude and frequency in order to fluidize the granular 
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bed to minimize resistance while keeping the bed rigid 
enough that it can support applied loads. For millimeter- 
sized undulatory swimmers, such as the worm nematode 
Caenorhabditis elegans {C. elegans)^ experiments show 
that the nematode swims more efficiently in a monodis- 
perse particulate system of packing fraction (j) = 0.6 than 
in a fluid without particles [13]. This observed swim- 
ming enhancement was found to be independent of the 
surrounding particle size. 

Despite recent efforts, many important questions re- 
main. For instance, it is known that certain properties 
of granular media, including packing fraction and particle 
size dispersity, can significantly affect the response of the 
bulk material to external perturbations due to the struc- 
ture of force chains. For example, the response of static 
quasi-2D piles to an external force behaves more rigid with 
increasing ^, and for disordered systems the response is 
best described by elastoplastic models [T5HT7] . It is yet 
to be seen how the difference in information propagation 
between ordered and disordered granular media can affect 
the motility of live organisms. 

Here, we investigate the motility behavior of the ne- 
matode C. elegans in a buffer solution ^ of viscosity 
and density similar to water. The solution contains ei- 
ther monodisperse or polydisperse particles of varying area 
density. Such a fluid can exhibit both solid- and fluid-like 
behavior depending on packing fraction cf) (referred to here 
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Fig. 1: (Color online) (a) Top view schematic of nematodes 
swimming in buffer solution (24 x 40 mm) and a quasi- 2D gran- 
ular layer (12 x 12 mm), (b) Side view schematic of the plate 
spacing used to confine the nematodes, buffer solution, and 
granular monolayer. The largest particles used were just under 
100 /im while the plate spacing was controlled by the volume 
of buffer solution and was approximately 300 /xm (c) Plot of 
the size distribution for monodisperse particles (out of 50 mea- 
surements) is shown in red, and for poly disperse particles (out 
of 100 measurements) is shown in gray. 



as particle area density) [19 and is believed to mimic more 
accurately the soil-like environments where the nematode 
is originally found. C. elegans is studied extensively and 
serves as a model organism, and are chosen for this study 
because they are known to move in different fluidic en- 
vironments [2QH23] . They are typically 1 mm in length, 
0.07 mm in diameter, and are equipped with 95 muscle 
cells aligning their dorsal and ventral sides [24] . Motility is 
achieved from a 'thrashing' motion against the surround- 
ing medium [25^, and typically occurs at a low Reynolds 
number. Re ~ O(10~^) [26]. This makes C. elegans an 
attractive candidate for motility studies. 

Experimental methods. — Experiments are per- 
formed in a shallow two-dimensional particle bed in order 
to minimize out-of-plane undulations. The granular mono- 
layer is created by placing the nematode, buffer solution, 
and particles between two glass cover-slips (170 /im thick), 
as shown in fig.[T]^a,b). The plate spacing is set by the con- 
tact angle between the liquid, the glass cover-slip, and the 
ambient air as well as the volume of the buffer solution. 
A constant volume of 400 (A of buffer solution is used for 



each experiment to produce a plate spacing of approxi- 
mately 300 /im, which is optically measured. The parti- 
cles then settle on the bottom plate creating a randomly 
packed granular layer in which the nematode moves free 
of confinement effects from the top cover-slip, as shown 
in fig. mb). Note that the total bath size (24 x 40 x 0.3 
mm) is much larger than the granular bed (12 x 12 x 0.08 
mm) and the nematode. Monodisperse samples consist 
of glass beads with diameter dm = 60 ± 3 /im (Mo-Sci 
Specialty Products) while polydisperse samples consist of 
glass beads with diameter dp = 52 ± 10 /im (Potters In- 
dustries Inc.). A histogram with the measured size distri- 
butions for both samples is shown in fig. Hj^c). 

Two main types of experiments are performed, body 
tracking methods and particle tracking velocimetry. The 
motility behavior of C. elegans is characterized using body 
posture tracking methods [27 . The effects of the nema- 
tode's motion on the surrounding particles are investigated 
using particle tracking methods [28 . Both C. elegans and 
particles are imaged using an epi-fluorescent inverted mi- 
croscope at 5x magnification (Axio Observer Zl, Carl 
Zeiss Inc.) with an acquisition rate of 125 frames per sec- 
ond (Fastcam SAl.l, Photron USA Inc.). The C. elegans 
strain (Genetic Stock Center) is cultured using standard 
methods [T8] . 

Nematode kinematics. — The C. elegans^ kinemat- 
ics in wet granular media is analyzed by tracking the ne- 
matode body posture in both space and time. Figure[2]and 
supplementary online material shows sample snapshots of 
the nematode moving in monodisperse (left) and poly- 
disperse (right) wet granular media and the correspond- 
ing body postures at a particle area density of approxi- 
mately 0.55. An important quantity that can be obtained 
from body postures is the undulation speed defined as 
A/, where A and / are the nematode's swimming ampli- 
tude and frequency, respectively. The beating amplitude 
A is defined as half of the peak-to-peak distance swept 
out by the nematode head perpendicular to the swimming 
direction (fig. [21 top). The beating frequency is the in- 
verse of the period T, or the time it takes for the head 
to complete one sweeping motion. The quantity Af is of- 
ten used as a relative measure of changes in motility gait 
[25 . Figure [3]^a) shows Af as a function of local area den- 
sity for both monodisperse and polydisperse media. We 
define the local area density (j) of the granular layer as 
(j) = N[d/3Ln]'^ ^ where N is the total number of spheres, 
with average diameter d of either 62 jam for monodisperse 
media or 52 jam for polydisperse media, within an interro- 
gation region locally around the nematode centroid with 
radius of 7^ = 3Ln/4, where Ln is a nematode body length. 

Results show that the undulation speed decreases by al- 
most one-half in wet granular media for area density rang- 
ing from 0.3 to 0.5. The undulation speeds are 0.52 ±0.05 
mm/s and 0.23 ± 0.05 mm/s for buffer solutions {(f) = 0) 
and granular media (mono- and polydisperse media), re- 
spectively. This decrease in Af indicates that the nema- 
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Fig. 2: (Color online) Plots of the nematode skeleton for 
two environments with particle area densities of 0.55; (Left) 
monodisperse granular media (Right) polydisperse granular 
media. The undulation period T is color-coded such that blue 
and red correspond to T = and T = 27r, respectively. Be- 
low the contours are time-sequenced micrographs of wild-type 
C. elegans moving through (Left) monodisperse and (Right) 
polydisperse wet granular media with an area density of 0.55. 
Images are taken at equally spaced intervals of At = 1.12 s pro- 
gressing from top to bottom. Note the differences between the 
bottom two micrographs: (Right) large amount of white space, 
or void space, free of particles compared to (Left) monodisperse 
media. See supplementary online material for corresponding 
movies. 

tode moves with slower and/or smaller undulations. For 
values of (j) larger than 0.55, the undulation speed is also 
affected by the media dispersity (fig. [3^). For the polydis- 
perse case, the undulation speed decreases for (j) > 0.55, 
whereas Af remains nearly constant for the monodisperse 
case under similar conditions. 

The decrease in undulation speed for area densities 
larger than 0.55 is accompanied by the nematode's change 
in motility gait. The change in motility behavior can be 
seen at the top of fig. [21 which shows the time-evolutions 
of nematode skeletons over one beating cycle. The contour 
lines are evenly spaced with intervals of 0.08 s, centered 
at the nematode's center-of-mass, with the head oriented 
to the right. The nematode exhibits markedly different 
waveforms in a buffer solution without particles (0 = 0) 
than in the same solution containing a high concentration 
of polydisperse particles {(j) > 0.55). For the monodis- 



Fig. 3: (Color online) (a) Undulation speed and (b) swimming 
speed as a function of local area density for C. elegans in wet 
granular media. The undulation speed decreases drastically in 
granular media for area densities above 0.3 and remains con- 
stant in monodisperse media, but decreases again in polydis- 
perse media beyond ~ 0.55. The swimming speed increases 
slightly in granular media and remains nearly constant with 
increasing 0. 



perse case, the nematode exhibits a waveform identical to 
the buffer solution {(j) = 0) for all values of (j) investigated 
here. The near-symmetrical waveform seen during swim- 
ming (fig. [21 top-left) in a buffer solution or monodisperse 
media is transformed into an apparent slithering motion 
as the nematode responds to a dense polydisperse medium 
(fig. [21 top-right). The peak-to-peak amplitudes for the 
two cases are approximately 0.36 mm for ^ = and 0.175 
mm for (j) > 0.55. This change in gait was not observed 
in monodisperse media, and it is indicative of the differ- 
ences in response and material properties between the two 
granular media, discussed in more detail in the following 
section. 

Figure [3](b) shows the nematode's propulsion speed U 
as a function of <p for monodisperse and polydisperse me- 
dia. Results show that U increases in wet granular media 
compared to swimming in a buffer solution, and remains 
relatively constant, even as (j) increases. The average ne- 
matode propulsion speed in buffer solution is 0.36 ± 0.01 
mm/s, whereas the propulsion speed in wet granular me- 
dia is 0.48 =b 0.12 mm/s, varying slightly with local area 
density. Overall, we find that even though the undulation 
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Fig. 4: (Color online) Stroke efficiency, the ratio of propulsion 
speed IJ to undulation speed A/, as a function of local par- 
ticle area density for C. elegans in wet granular media. The 
Stroke efficiency increases with area density up to 0.55 and then 
remains constant at a value of approximately 2.25 for monodis- 
perse media. The Stroke efficiency increases monotonically for 
polydisperse media. 



speed decreased in dense wet granular media, this did not 
hinder the C. elegans ability to move. This observation 
can be further quantified by computing the ratio of propul- 
sion speed U to undulation speed Af as a function of (/), 
referred to here as the Stroke efficiency u = U/Af^ shown 
in fig. m Values of u > 1 imply that there is greater for- 
ward swimming speed than transverse undulation speed. 
We find that u ~ 0.69 for swimming in a fluid without par- 
ticles and u ~ 2.08 for swimming in wet granular media. 
The Stroke efficiency increases with increasing cj) implying 
that the nematode moves more efficiently by generating 
more forward motion with less lateral motion. Obviously, 
the efficiency cannot be unbounded and is expected to 
plateau at some value of (j) until reaching the jamming 
limit. Monodisperse media appears to plateau to a value 
of z/ ^ 2.25 for area densities greater than 0.55. We note 
that the Strouhal number, which is the inverse of is 
typically used to characterize vortex shedding in oscillat- 
ing flow mechanisms. It was also recently used to suggest 
that C. elegans move more efficiently in granular media 

Nematode propulsion in granular media. — It is 

known that C. elegans have distinct motility gaits, such 
as swimming or crawling, depending on whether the en- 
vironment is a liquid buffer solution or a solid agar gel, 
respectively [25 . That is, the nematode is sensitive to the 
material properties of its surrounding fluidic environment. 
Here, for (j) < 0.3, particles are sparse (dilute) and the ne- 
matode behaves essentially as if it swimming in a pure 
buffer solution with no particles. In addition, figs. [3] and |4] 




Fig. 5: (Color online) Particle pathlines induced by the ne- 
matode in (a) monodisperse and (b) polydisperse media for 
(/) 0.55. The corresponding nematode undulation frequencies 
are fm ~ 1.5 Hz and /p 0.5 Hz. Particles are tracked for 
t = 3.1 s. Nematode trajectories are depicted by the blue line 
with circles representing the nematode center of mass. The 
length of particle trajectories is color coded such that white 
and red correspond to short and long trajectories, respectively. 
The difference in rearrangement to nematode undulations is ev- 
ident. Monodisperse particle trajectories are localized around 
the main nematode trajectory, whereas polydisperse particle 
trajectories extend outward normal to the main nematode tra- 
jectory. See supplementary online material for corresponding 



show that for area densities in the range of 0.3 < (\) < 0.55, 
the nematode motility behavior is similar for both mono- 
and polydisperse media which indicates that both media 
have similar material properties. This is due to the lack 
of large scale structure formation in both media for this 
semi-dilute density range. However, for (j) > 0.55, there is 
an obvious difference in motility behavior that is specific 
to mono- or polydisperse media as noted by the wave form 
(fig. [21 top) and undulation speed (fig. [3|). This suggests 
that in this dense regime (^ > 0.55) the bulk material 
properties of polydisperse media are similar to an elastic 
material [1 6,29.30 such as an agar gel, while monodisperse 
media behave similar to a viscous fluid [31 . 

In order to gain further insight into the difference in re- 
sponse of mono- and polydisperse media to the nematodes 
undulations for (j) > 0.55, we tracked an average of 2,000 
particles along the nematodes path. Figure [5] shows par- 
ticle pathlines for monodisperse (fig. [5^) and polydisperse 
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Fig. 6: (Color online) (a) Average overall particle displacement as a function of the radial distance away from the nematode body; 
the gray region corresponds to the body thickness. Particle displacement is largest near the nematode body and decays rapidly 
for monodisperse media compared to the polydisperse media, (b) PDF of particle velocities localized around the nematode 
trajectory. While particle trajectories are longer for monodisperse media (see fig.Ej), the probability density function (PDF) of 
the particle velocities are similar for both monodisperse and polydisperse media, (c) PDF of (mod tt), the angle between the 
nematode trajectory and surrounding instantaneous particle velocities. The difference in rearrangement between the two media 
is evident; the majority of monodisperse particles move nearly anti-parallel to the nematode main trajectory while polydisperse 
particles move perpendicular to the nematode main trajectory. Note that the lines are added to guide the eye. 



(fig.ISb) media as a nematode moves through a 3.1 mm by 
1.8 mm region for (j) ^ 0.55. For the monodisperse case, 
most of the long particle trajectories are localized around 
and are tangential to the nematode's body. By contrast, 
for the polydisperse case, particle displacements extend 
further away from and are mostly normal to the nema- 
tode's body. The average displacement of particles near- 
est the nematode is Ax^ ^ 60 /im and Axp ~ 40 /im for 
monodisperse and polydisperse media, respectively. The 
average particle displacement decreases rapidly with radial 
distance from the nematode's body (fig. [6^); this decay is 
expected for viscous fluids. The average particle displace- 
ment decays to one-half of its magnitude at a distance of 
0.3Ln for monodisperse media compared to a distance of 
0.55Ln for polydisperse media, where = 1 mm. As 
expected, the distributions of particle velocities for both 
media are similar, which indicates similar energy input 
from nematode undulations (fig. [GJd). 

In general, we find that there is a noticeable difference in 
particle rearrangement behavior between mono- and poly- 
disperse media due to nematode undulations at compara- 
ble (j). We characterize this difference by computing the 
angle between the nematode trajectory and the instanta- 
neous velocity of the surrounding particles, A6> (mod tt), 
shown in fig. [6fc). Here, the values of A6> (mod tt) = 0° 
and A^ (mod tt) = 180° correspond to the same and oppo- 
site direction of swimming, respectively. The data shows 
that for the monodisperse medium, there is a sharp peak 
at 0.8 or A^ ^ 150° which indicates that the particles 
move along the nematode's main trajectory. By contrast, 
for the polydisperse medium, there is a broad distribu- 
tion of particle directions mostly moving normal to the 
nematode's main trajectory. The sample snapshots of the 
nematode moving in dense granular media (fig. [2] and sup- 



plementary online movies) illustrates the difference in re- 
arrangement behavior. For the monodisperse case, we ob- 
serve that the local area density at the tail is only slightly 
less than at the head (fig. [21 bottom-left) indicating that 
particles are able to easily rearrange during the perturba- 
tion imposed by the moving nematode. This behavior is 
similar to what is observed for an intruder dragged at con- 
stant force through an amorphous monolayer of vibrated 
grains [32 . For the polydisperse case, however, the overall 
motion of both small and large particles is mostly laterally 
from the trajectory of a nematode (fig. [21 bottom-right) 
and particles do not readily fill the voids, leaving a dilute 
path behind the nematode. This observation provides fur- 
ther evidence that polydisperse media behaves elastoplas- 
tically, by retaining deformations due to nematode undu- 
lations. 

Overall, we find that C. elegans^ are able to adapt and 
utilize the surrounding granular media to convert lateral 
undulations into net forward motion. This conversion 
leads to an increase in propulsion speed, which the ne- 
matode maintains even with increasing particle area den- 
sity and degree of size dispersity. The observed enhanced 
propulsion may be explained by closely inspecting the 
nematode skeleton during swimming in buffer solutions, 
(fig. [21 top- left), and recognizing that the motion over 
one beating cycle is mostly in the lateral direction. As 
(j) increases, the surrounding particles act as supporting 
structures that convert undulatory motion into net for- 
ward motion. Here, the effective propulsion depends on 
the difference between the applied load by the nematode 
undulations and the rigidity of the granular media ^ . 
The rigidity, or ability to withstand deformation, depends 
on the local particle area density and becomes more solid- 
like at larger densities, until reaching the jamming limit 
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~ 0.84 [33^. This rigidity also depends on the size dis- 
persity, which affects the amount of disorder of the gran- 
ular media [T7l[34]. Polydisperse samples are more disor- 
dered and therefore present a more "rigid" quasi-static flu- 
idic environment, whereas monodisperse samples are less 
disordered and are more easily sheared, therefore appear 
more fluid-like to the C. elegan [35 . 

Finally, an estimate of the propulsive force {Fp) can be 
obtained by computing the work done on the surrounding 
particles from nematode undulations, W = FpA — Fdx, 
where A is the nematode beating amplitude, Fd is Stokes' 
drag, and x is the average displacement of a group of par- 
ticles. The kinetic energy of the particles is computed 
from particle tracking data. The propulsive force is es- 
timated by Fp = 7rpd^[Nv'^g/12A] + 37rr]d[Nvgx/A]. Here 
p is the particle density, d is the particle diameter, r] is 
the buffer dynamic viscosity, N is the average number of 
particles moving, and Vg is the average group velocity of 
the particles per undulation. This approximation reveals 
that the force of the C. elegans on the surrounding media 
is Fp ~ (9(10~^^) Newtons. This value is similar to the 
force generated by a slender-body swimming in a viscous 
fluid [26] and is only slightly less than values measured 
for crawling C. elegans on substrates using force sensing 
pillars [36 . 

Conclusions. — In conclusion, we have shown that 
live organisms can exploit the material properties of the 
surrounding complex media to achieve efficient propulsion. 
Surprisingly, the addition of particles to a buffer solution 
does not hinder the C. elegans ability to move. In fact, the 
nematode's propulsion speed is enhanced by the presence 
of particles and is not dependent on particle area density. 
The size dispersity on the other hand, does affect the lo- 
comotion gait of the nematode at high area densities. 

The adaptive motility behavior of C. elegans in wet 
granular media may explain the consistent performance 
seen not only in nematodes [20] , but also in flagellated bac- 
teria and human spermatozoa. Similar mechanisms have 
been suggested in the enhanced movement of Leptospira 
[37] and other flagellated bacteria in non-Newtonian flu- 
ids. Hence, it is possible that other forms of microscopic 
propulsion are also tuned to biological fluids which pos- 
sess anisotropics. In addition, we are able to study the 
response and rearrangement of the granular media and its 
dependence on packing fraction and particle size disper- 
sity free of external forces, such as gravity. This suggests 
an interesting method to probe more complex anisotropic 
media that could depend on shape or particle interactions. 
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